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ABSTRACT 

We present the results of a A = 2.7 mm continuum interferometric survey of 24 
young stellar objects in 11 fields. The target objects range from deeply embedded 
Class sources to optical T Tauri sources. This is the first sub-arcsecond survey of 
the A = 2.7 mm dust continuum emission from young, embedded stellar systems. 
These multi-array observations, utilizing the high dynamic u,v range of the BIMA 
array, fully sample spatial scales ranging from 0'.'4 to 60", allowing the first consistent 
comparison of dust emission structures in a variety of systems. The images show 
a diversity of structure and complexity. The optically visible T Tauri stars (DC 
Tauri, HL Tauri, GG Tauri, and CM Aurigae) have continuum emission dominated 
by compact (< 1") circumstellar disks. In the cases of HL Tauri and DG Tauri, 
the disks are resolved. The more embedded near-infrared sources (SVS13 and L1551 
IRS5) have continuum emission that is extended and compact. The embedded sources 
(L1448 IRS3, NGC1333 IRAS2, NGC1333 IRAS4, VLA 1623, and IRAS 16293-2422) 
have continuum emission dominated by the extended envelope, typically > 85% of the 
emission at A = 2.7 mm. In fact, in many of the deeply embedded systems it is difficult 
to uniquely isolate the disk emission component from the envelope extending inward 
to AU size scales. Simple estimates of the circumstellar mass in the optical/infrared 
and embedded systems range from 0.01-0.08 MQand 0.04-2.88 Mq, respectively. All 
of the target embedded objects are in multiple systems with separations on scales 
of ~ 30" or less. Based on the system separation, we place the objects into three 
categories: separate envelope (separation > 6500 AU), common envelope (separation 
150-3000 AU), and common disk (separation < 100 AU). These three groups can be 
linked with fragmentation events during the star formation process: separate envelopes 
from prompt initial fragmentation and the separate collapse of a loosely condensed 
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cloud, common envelopes from fragmentation of a moderately centrally condensed 
spherical system, and common disk from fragmentation of a high angular momentum 
circumstellar disk. 

Subject headings: stars: formation — stars: pre-main-sequence — survey — binaries: 
general — radio continuum: stars — infrared: stars — techniques: interferometric 

1. Introduction 

Young stellar systems exhibit excess infrared and millimeter emission that arises primarily 
from circumstellar dust in two basic evolutionary structures: envelopes and disks. Current 
observations and theories of star formation (e.g. Larson 1969; Penston 1969; Shu 1977; Cassen 
k Moosman 1981; Lada & Wilking 1984; Adams, Lada, & Shu 1987; Shu et al. 1993; Andre, 
Ward- Thompson, & Barsony 1993) outline an evolutionary sequence that begins with a density 
enhancement which quasi-statically contracts to form a centrally concentrated core. The core then 
gravitationally collapses forming an infall region which feeds a central protostar and a surrounding 
disk supported by centrifugal forces. Initially, the envelope (radii of 1000s of AU) contains most 
of the mass, but as the system evolves, the disk and protostar grow, and the disk becomes a 
significant mass reservoir. In time, the prenatal envelope is depleted of material and progressively 
blown away by the powerful stellar outflow, revealing the young star and disk (radii of ~ 100 AU) 
system. Detailed modeling of young stellar objects has shown that the infrared through millimeter 
wavelength spectral energy distribution (BED) can be reproduced by thermal dust emission from 
a combination of a large-scale envelope, spatially thin disk, and central star (Adams, Lada, &: Shu 
1987; Kenyon &; Hartmann 1987; Bertout, Basri, &; Bouvier 1988; Calvet, Hartmann, & Whitney 
1994) or, in some cases, by just a disk and central star (Beckwith et al. 1990; Osterloh & Beckwith 
1995; Dutrey et al. 1996). 

Surveys of main-sequence stellar systems find that most stars are in binary or multiple 
systems (Duquennoy &; Mayor 1991) with the separation distribution ranging from a few Rq to 
10^ AU and the probability distribution peaking around 30 AU. Recent surveys of the nearby 

star-forming regions of Taurus and Ophiuchus find that the occurrence of binaries in the young 
visible systems is about twice as common as among main-sequence stars (Simon et al. 1992; Ghez, 
Neugebauer, &, Matthews 1993; Leinert et al. 1993; Reipurth &: Zinnecker 1993; Ghez, White, 
&; Simon 1997). How is the above star formation process altered to form binary systems, and 
how does the younger, deeply embedded systems binary occurrence compare to the young visible 
systems? 

In this paper, we present a A = 2.7 mm continuum survey of 24 nearby Young Stellar Objects 
that represent a sample of young stellar systems at various stages of evolution. The survey 
highlights the large dynamic range of u,v spacings available with the Berkeley-Illinois-Maryland 
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Association (BIMA) millimeter array — covering both the largest and smallest u,v spacings 
currently available at A = 2.7 mm. With a combination of low and high rcsohitions, we are able 
to map the envelopes of the embedded sources and resolve-out the large-scale structure in order 
to peer inside the envelopes and image the central regions. With the highest angular resolution 
to date at this wavelength, we are able to image circumstellar disks (e.g. Mundy et al. 1996) 
and search for close binaries (e.g. Looney, Mundy, & Welch 1997). The purpose of this paper 
is to present our images with the discussion focusing on differences and similarities between the 
various evolutionary stages and several broad trends in the data. In additional papers, we will 
discuss individual sources in detail and extensively model the emission structures of these sources 
as arising from extended envelopes, circumstellar disks, and circumbinary disks. 

2. Sample, Observations, and Mapping 
2.1. Sample 

The goal of the survey was to image a range of young stellar systems at sub-arcsecond 
resolution. We concentrated on bright, nearby sources which were most likely to produce high 
dynamic range images. The criteria for selecting the sample were: (1) for the best sensitivity 
to solar-system-sized spatial scales, we focused on the closest sources (< 350 pc); (2) to insure 
sufficient surface brightness at sub-arcsecond resolution, we chose among the brightest sources 
in the nearby star forming regions; (3) to survey a number of evolutionary stages, our sample 
included the youngest, most deeply embedded sources (Class 0), as well as optical T Tauri stars 
(Class I/II); (4) finally, to achieve the best map fidelity, we required that the sources have both 
a bright phase reference calibrator and a weaker point source response calibrator nearby on the 
sky. This combination of criteria introduces three strong biases in our sample: (a) we favor bright 
millimeter sources which should be systems with more circumstellar material than average (e.g. 
Beckwith et al. 1990), (b) most of our sources are either Class or Class II (this is mostly due 
to the bias above), and (c) as a result of our distance limit, our sources are drawn from just 
three local clouds — Taurus, Ophiuchus, and Perseus. Table 1 lists the sources, adopted distances, 
whether they are optically visible (in this category we also include objects which are visible in 
the near-infrared) or embedded, adopted SED Class, main calibrator, secondary calibrator, and 
distance reference. The distance to the Perseus objects has been in dispute lately, 350 pc (Herbig 
Sz Jones 1983) and more recently 220 pc (Cernis 1990). We adopt the previous value, allowing an 
easier comparison to earlier works. 

The frequency of multiple systems will be addressed in §6, but the number of objects in the 
embedded source group is misleading since all of the observed embedded sources are in multiple 
systems. For example, the Perseus objects include 13 sources which were mapped in only 4 
pointings, whereas in the optical group, only DG Tauri and DG Tauri B (which is an embedded 
object not thought to be related to DG Tauri except by projected proximity; Jones & Cohen 1986) 
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were mapped in the same observation. 

2.2. Observations 

All sources were observed in three configurations (C, B, and A) of the 9-element BIMA 
Array|(Welch et al. 1996). The observations were acquired from 1996 May to 1998 March. The 
digital correlator was configured with two 700 MHz bands centered at 106.04 GHz and 109.45 GHz. 
The 0^*^0(1-0) line was observed simultaneously; those results will be discussed elsewhere. The 
two continuum bands were checked for consistency, then combined in the final images. The system 
temperatures during the observations ranged from 150-700 K (SSB). 

In the compact C array (typical synthesized beam of ~8"), the shortest baselines were limited 
by the antenna size of 6.1 m, yielding a minimum projected baseline of 2.1 kA and good sensitivity 
to structures as large as ~50". In the mid-sized B array (typical synthesized beam of ~2"), the 
observations are sensitive to structures as large as ~10". In the long baseline A array (typical 
synthesized beam of ~0'.'6), the longest baselines were typically 450 kA. The combination of these 
three arrays provides a well sampled u,v plane from 2.1 kA out to 400 kA. 

The uncertainty in the amplitude calibration is estimated to be 20%. In the B and C arrays, 
the amplitude calibration was boot-strapped from Mars. In the A Array, amplitude calibration 
was done by assuming the flux of the quasar 3C273 to be 18.8 Jy at the end of 1996 October, 23.0 
Jy at the end of 1996 November, and 27.0 Jy until the end of 1997 January. Absolute positions 
in our map have uncertainty due to the uncertainty in the antenna locations and the statistical 
variation from the signal-to-noise of the observation. These two factors add in quadrature to give 
a typical absolute positional uncertainty of O'.'IS in the highest resolution maps. 

The A array observations required careful phase calibration. On long baselines, the 
interferometer phase is very sensitive to atmospheric fluctuations. We employed rapid phase 
referencing; the observations were switched between source and phase calibrator on a two minute 
cycle, to follow the atmospheric phase (Holdaway & Owen 1995; Looney, Mundy, & Welch 1997). 
To monitor the coherence of the atmosphere, or the "seeing" , we included a nearby weaker quasar 
in the observation cycle. This quasar was imaged along with the target source as a check of the 
point source response and for accurate image registration. In the observations presented here, the 
secondary quasar always mapped as a point source within statistical uncertainties. 



The BIMA Array is operated by the Berkeley lUinois Maryland Association under funding from the National 
Science Foundation. 
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2.3. Displaying the Data 

The observational data span u,v distances from 2.1kA to 450kA, providing information of the 
brightness distribution on spatial scales from (X'4 to 60". In order to display this information in 
the image plane, we have mapped the emission with four different u,v weighting schemes which 
stress structures on spatial scales of roughly 5", 3", 1", and (/.'G. These resolutions were typically 
obtained with natural weighting, robust weighting (Briggs 1995) of 1.0, robust weighting of 0.0, 
and robust weighting of -0.5, respectively. 

3. Results 

The A = 2.7 mm continuum images from the survey are shown in Figures 1 through 16, first 
the optical/IR systems and then the deeply embedded systems. In each figure, the four panels 
display the same multi-configuration data with different u,v weighting schemes to emphasize 
structures on size scales comparable to the synthesized beam. Table 2 lists, at each resolution, the 
peak flux, the integrated flux, and the box used to measure the integrated flux. The error bars on 
the flux measurements represent the statistical uncertainties. 

In Figure 17, the source fluxes are presented in a plot comparing the total integrated flux 
to the ratio of the visibility amplitude at two specific fringe spacings. The horizontal axis is the 
integrated flux of each object taken from Table 2 and adjusted to the distance of Taurus (140 
pc). For the vertical axis, the u,v data were binned in annuli stretching from 2.2 kA to 7.8 kA 
and from 22 kA to 78 kA, corresponding to average spacings of 5 kA and 50 kA for the distance of 
Taurus. The ranges were chosen to provide the best signal to noise for a typical source and the 
best sensitivity to large and small scale structures. The vertical axis plots the ratio of the vector 
averaged amplitudes in the two bins or the 5kA/50kA ratio. In order for the ratio to consistently 
probe the same spatial scales for all objects, we adjusted the bin range to take into account the 
source distance: for p Ophiuchi we used annuli averaging of 5.7 kA and 57 kA, respectively, and 
for Perseus we used annuli averaging 12.5 kA and 125 kA, respectively. For objects in multiple 
systems, the companion objects were subtracted out of the u,v data using the CLEAN components 
in the lower resolution (5" and 3") maps as a model. 

Figure 18 illustrates how to interpret Figure 17 by comparing the visibility curves of a 5000 
AU outer radius envelope and a 100 AU radius disk, both at the distance of Taurus (140 pc). 
This figure is essentially the one-dimensional Fourier Transform of the brightness distribution of 
the source on the sky. The horizontal axis is the spatial frequency of the interferometer, or the 
antenna separation in units of wavelengths, and the vertical axis is in Janskys. The visibility of 
the disk model is similar to an unresolved point source until the disk starts to be resolved at ~ 
60 kA. In the case of the envelope, it is quickly resolved and approaches the expected power-law 
form; the "wiggles" are due to the sharp outer edge of the model. Overall, the 5kA/50kA ratio 
provides a measure of the relative emission on spatial scales of 2500 - 9000 AU and 250 - 900 AU 
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and quantitatively measures the "compactness" of the object. An object with spatially extended 

structure, such as an envelope with size scales of 1000 AU or larger will have a 5kA/50kA ratio >1, 
and an object that is entirely compact, such as a circumstellar disk with radius of ~ 100 AU, will 
be essentially unresolved in the two bins and the ratio will be unity. 

The general trend of Figure 17 is, as expected, that most of the optical sources (solid 
triangles) are compact (5kA/50kA ratio of ~ 1) and most of the embedded sources (solid squares) 
are surrounded by envelopes that are being significantly resolved at 50 kA (5kA/50kA ratio >1). 
However, there are a couple of exceptions worth discussing. First, there are two optical stars 
with unusually large 5kA/50kA ratios — GG Tauri and SVS13 A. GG Tauri is a close binary 
system with a separation of 0'.'255 (Leinert et al. 1991) and a large circumbinary disk (diameter 

400 AU; Simon & Guilloteau 1992; Dutrey, Guilloteau, & Simon 1994). Thus, in GG Tauri 
the 5kA/50kA ratio is resolving the large scale circumbinary disk. SVS13 A, first detected in the 
infrared at 2.2 fim (Strom, Grasdalen, &; Strom 1974; Strom, Vrba, & Strom 1976) is also know 
to have optical/infrared outbursts (Eisloffel 1991); yet it has a very large 5kA/50kA ratio. There 
are two possible explanations: (1) the envelope of the nearby, younger embedded object SVS13 
A2 is contributing to the flux of SVS13 Al at 5 kA, or (2) the optical/infrared emission is a 
reflection nebula and the source should be classiflcd as embedded. A second set of exceptions are 
two embedded sources (NGC 1333 IRAS2-B and IRAS4-C) that have unusually smah 5kA/50kA 
ratios. Since these sources are not detected in the optical or the near-infrared they are classified 
as embedded sources, but their 5kA/50kA ratio and their integrated fluxes in Table 2 indicate 
that they arc compact. These two sources could be optical/near-infrared sources that are viewed 
through intervening obscuration. The following subsections discuss each of the sources in more 
detail. 



3.1. DG Tauri 

DG Tauri is a well studied classical T Tauri star system. Through modeling of the system's 
SED, Adams, Emerson, & Fuller (1990) estimated a radius of 75 AU for the circumstellar disk. 
The source was observed in the near-infrared during a lunar occultation (Leinert et al. 1991), and 
it was determined that the star was a single system with an extended "shell" 6.8 AU in diameter. 
In addition, near-infrared speckle observations revealed the presence of a "halo" with a diameter 
of 130 AU (Leinert et al. 1991). In panel (d) of Figure 1, the circumstellar disk around DG 
Tauri is marginally resolved. Fitting an elliptical Gaussian to the image in panel (d). we obtain 
a deconvolved Gaussian size of O'.'Gl it O'.'l x 0'.'57 it O'.'l with a principal axis of 167° ± 10°. 
This result is different from the measurement at A = 2.0 mm from Kitamura, Kawabe, & Saito 
(1996), which found a deconvolved size of 1''56 x 0'.'54 at a principal axis of 99°, but similar to 
the measurement at A = 2.7 mm from Dutrey et al. (1996), which found a deconvolved size of 
I'.'l X 0'.'6 at a principal axis of 150°. The extension to the southwest in panel (d) lies along the 
direction of the jet (Cohen & Bieging 1986; Stapelfeldt et al. 1997; Lavalley et al. 1997) and may 



- 7- 



trace ionized gas in the jet. 

3.2. DG Tauri B 

DG Tauri B, located 53" southwest of DG Tauri, has a molecular outflow (principal axis 
of ~ 295°; Mitchell et al. 1994; Mitchell, Sargent, & Mannings 1997) that is driven by a jet 
seen at optical (Mundt, Brugel, & Biihrke 1987) and centimeter (Rodriguez, Anglada, & Raga 
1995) wavelengths. DG Tauri B was observed near the half power point of our beam during the 
observation of DG Tauri; therefore measured fluxes have a significant additional uncertainty. The 
fluxes listed in Table 2 were corrected for the primary beam attenuation. 

The morphology of DG Tauri B changes in Figure 2 with increasing resolution. Going from 
panel (d) to panel (c) to panel (b), the major elongation of the emission changes from northwest to 
north to slightly northeast. In panel (a) the emission is triangular with extension to the northwest, 
northeast, and southwest. The simplest interpretation is that the high resolution image traces the 
ionized gas, while the lower resolution images trace both ionized gas and dust. The position angle 
for the larger scale dust emission is then ~35°, which is consistent with the optical extinction lane 
(Stapelfeldt et al. 1997) and perpendicular to the outflow jet. The extended emission in panel (d) 
strongly resembles the A = 3.6 cm image (Rodriguez, Anglada, h Raga 1995), suggesting that it 
is tracing ionized gas in the jet. The relative flux numbers in Table 2 suggest that roughly half of 
the flux arises from dust and half from ionized gas in the jet. 

3.3. L1551 IRS5 

L1551 IRS5, one of the prototypical Class I sources in the classification scheme of Adams, 
Lada, & Shu (1987), has the most spectacular bipolar molecular outflow in the Taurus cloud 
(principal axis of ~ 50°; Snell, Loren, & Plambcck 1980). The A = 2.7 mm continuum data 
presented here were discussed in detail by Looney, Mundy, & Welch (1997). The source is argued 
to be a proto-binary system with a large-scale envelope (~1300 AU radius), circumbinary disk, 
and two circumstellar disks (separation of 0'.'35). The binary circumstellar disks have recently been 
resolved in A = 7 mm observations (Rodriguez et al. 1998). In Figure 3 panels (a) and (b), the 
emission is dominated by the large-scale envelope, while panel (c) clearly shows the circumbinary 
envelope. In panel (d), the two point-source-like circumstellar disks are still convolved with the 
low-level emission from the circumbinary envelope which is extended along a principal axis of 
~ 160°. The higher resolution image from Looney, Mundy, & Welch (1997) is not shown. 
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3.4. HL Tauri 

HL Tauri, perhaps the most studied of the optical/IR visible young stars, has a large-scale 
CO structure 1500 AU, Sargent & Beckwith 1991; Hayashi, Ohashi, & Miyama 1994) and a 
compact circumstellar disk (~ 100 AU) that has been resolved by the CSO-JCMT interferometer 
(Lay ct al. 1994; Lay et al. 1997) and imaged by the BIMA array (Mundy ct al. 1996). Figure 4 
shows the new BIMA image which has both lower noise and higher resolution than the image of 
Mundy et al. (1996). 

In panel (d), the circumstellar disk of HL Tauri is clearly evident. Fitting an elliptical 
Gaussian to the image, we obtain a deconvolved Gaussian size of O'.'SSztO'.'l x 0'.'58±0'.'1 and 
principal axis of 111°±10°, which agrees with the observations of Lay et al. (1994) and Mundy 
et al. (1996). However, fitting an elliptical Gaussian to the image is not the correct method 
for determining the true disk size. Recent modeling of the HL Tauri circumstellar disk found 
that simple models could not fit the CSO-JCMT single baseline interferometer A = 650 /xm and 
850 fim data and the A = 2.7 mm and 7 mm data (Lay ct al. 1997). We will consider more 
complicated disk models in a subsequent paper. The image in panel (d) also shows an extension 
to the north-east along the axis of the optical jet, principal axis 46° (Mundt et al. 1990). This 
feature likely arises from free-free emission in the jet; such free-free emission dominates the highest 
resolution maps at A = 7 mm (Wilner et al. 1999). 

HL Tauri is classified as an optical source, but has recently been shown to be embedded 
within a reflection nebula (Stapelfeldt et al. 1995); we do not see the star directly in optical light, 
but it can be seen directly in the near-infraxed (Weintraub, Kastner, Sz Whitney 1995; Beckwith 
& Birk 1995). Our data do not conclusively detect envelope emission associated with the extended 
nebula. The envelope on size scales larger than 3" contributes less than 10% of the dust emission, 
where our estimate is limited by the uncertainty of the relative amplitude calibration between 
arrays . 

3.5. GG Tauri 

GG Tauri is a close binary system with a separation of 0'.'25 (Leinert et al. 1991) and a large 
circumbinary "ring-like" disk (diameter ~ 400 AU; Simon &; Guilloteau 1992; Dutrey, Guilloteau, 
& Simon 1994). Our images presented in Figure 5, have different u,v weighting schemes from 
the rest of the surveyed objects stressing size scales of 5", 2", 1", and O'.'O. Fitting an elliptical 
Gaussian to the image in panel (b), we obtain a deconvolved size of 3'.'3 it O'.'l x 2'.'7 ± O'.'l at a 
position angle of 82° it 10°, which is in good agreement with Dutrey, Guilloteau, &: Simon (1994). 

Most of the peaks and valleys in the emission in panel (d) represent <2.5cr variations from an 
overall flat structure; hence, it is difficult to say whether they are real variations in the surface 
density or temperature. The emission is not strongly concentrated toward the inner edge of 
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the ring indicating that there is not a strong surface density gradient within the disk. A more 
detailed discussion of the disk will be presented in Mundy, Looney, &; Welch (1999). No emission 
is detected from the system at 0'.'6 resolution; this places upper limits on the emission from any 
compact structures (< C/.'G), such as individual circumstellar disks within the binary system, at a 
3(7 limit of 5 mJy. The companion binary system of this quadruple system, GG Tauri/c, was not 
detected at any resolution; the 3a upper limit on its flux density is 4 mJy. 

3.6. GM Aurigae 

GM Aurigae is another classical T Tauri star system that has a large-scale CO structure 
(Koerner, Sargent, & Beckwith 1993). In Figure 6 panel (d), we do not see evidence that the disk 
is resolved, but the signal-to-noise is poor. Fitting an elliptical Gaussian to the image in panel (d) 
yields a point source. We can place a limit on the deconvolved Gaussian size for the circumstellar 
disk of <0'.'4, at the 95% confidence level. In panel (c), the emission seems slightly extended along 
the direction perpendicular to the larger scale CO structure which has a position angle of 55°. 

The total flux density reported in Table 2 (22 mJy) is roughly consistent with that measured 
by Dutrey et al. (1996), 28 mJy Unlike Dutrey et al., we do not directly resolve the disk. 
However, we do measure a 35% decrease in flux density between the 5" and ff.'G beams, indicating 
that some structure is present. 

3.7. L1448 IRS3 Region 

The L1448 complex is located about ~1° southwest of NGC 1333. IRAS revealed three strong 
infrared sources in the region, of which L1448 IRS3 was the brightest in the far infrared (Bachiller 
& Cernicharo 1986). IRS3 is projected within the blueshifted lobe of the impressive, highly 
collimated molecular outflow from L1448-mm which lies to the southeast (Bachiller et al. 1990; 
Bachiller, Andre, &; Cabrit 1991). Coinciding within the uncertainties of the L1448 IRS3 source is 
a very strong H2O maser and A = 6 cm compact emission (Anglada et al. 1989). Higher resolution 
maps in the A = 2 cm and 6 cm continuum found that the source was composed of two sources 
L1448 N(A) and L1448 N(B) (Curiel et al. 1990; Barsony et al. 1998). Curiel et al. separated the 
region into two areas: L1448 C, the center of the molecular outflow, and L1448 N corresponding 
to the IRS3 source. L1448 N(B) contributes most of the flux at millimeter wavelengths (Terebey, 
Chandler, Sz Andre 1993; Terebey & Padgett 1997). A third source is also seen at A = 2.7 mm 
which lies to the north-west of L1448 N(B) (Terebey & Padgett 1997). 

In our images of the region, we clearly detect all three sources which we label: L1448 IRS3 A, 
B, and C, using the lAU nomenclature. The three sources are indicated in Figure 7 (b). Source 
A, which is the brightest source at centimeter wavelengths, is significantly weaker than source B 
at A = 2.7 mm. In fact at the highest resolution, source A is not detected. However, the spectral 
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index derived by Curiel et al. for source A (a ^ 0.2) suggests only ~ 2 mJy of the emission arises 
from free-free emission. Thus, the mihimeter wavelength flux from IRS3 A is likely dominated by 
dust emission. 

Located to the north-west, source C is detected at all resolutions. Unfortunately, source C 
is too weak to be plotted on the 5kA/50kA ratio figure. In panels (c) and (d), source B shows 
very complicated morphology on small scales. There is an outflow associated with the IRS3 region 
which is nearly parallel to the outflow from L1448-mm, at a position angle of ~ -21° (Bachiller 
et al. 1995; Davis & Smith 1995). The extension that is seen in panel (c) and (d) is almost 
perpendicular with the outflow, but it is unclear if it is an envelope or a large disk. The peak flux 
density decreases by a factor of two in each step of resolution in Figure 7(b), to (c), to (d), so 
there is considerable structure on all spatial scales. 

3.8. NGC 1333 IRAS2 

The NGC 1333 star forming region in Perseus is an extremely active site of star formation 
with multiple infrared sources (Strom, Vrba, Strom 1976; Aspin, Sandell, &: Russell 1994; Lada, 
Alves, & Lada 1996) and outflows (Sandell et al. 1994; Warin et al. 1996; Bally et al. 1996). NGC 
1333 IRAS2 (Jennings et al. 1987) is located on the edge of the large cavity in NGC1333 (Langer, 
Castets, & Lefloch 1996). The region has a two outflows that originate near IRAS2: the "N-S" 
outflow with principal axis of ~ 25° (Liseau, Sandell, & Knee 1988) and the "E-W" outflow with 
principal axis of ~ 104° (Sandell et al. 1994). Recent millimeter interferometric observations show 
that there are two continuum peaks that are probably associated with the two outflows, and that 
the northern source (Source A) is responsible for the "E-W" outflow (Blake 1997). 

Figures 8 and 9 show NGC 1333 IRAS2 A and B respectively. Source A, the stronger of the 
two sources, is mostly extended emission, but the flux remaining in panel (d) is consistent with 
a point source. Source B is mostly compact emission. The extension of source B in panel (d) is 
nearly perpendicular with the "N-S" outflow, suggesting a possible circumstellar structure. 

3.9. SVS 13 

Located northwest of IRAS2, the young stellar object SVS13 (Strom, Vrba, Strom 1976; also 
referred to as SSV13 in the literature from Herbig &; Jones 1983) is associated with the NGC 1333 
IRAS3 region (Jennings et al. 1987). IRAS3 is comprised of at least 3 millimeter sources: source 
A located near the infrared position of SVS13, source B to the southwest (Grossman et al. 1987; 
Chini et al. 1997; Bachiller et al. 1998), and source C further to the southwest (Chini et al. 1997; 
Bachiller et al. 1998). 

Figures 10 and 11 clearly show all three millimeter sources. In panel (b) of both figures 
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there is another source located to the southwest of source A. This source (which we will call 
A2) is coincident with VLA source 3 from recent VLA observations of this region (Rodriguez, 
Anglada, & Curiel 1997). Located ~ 6" from SVS13, Rodriguez et al. argue that A2 is a better 
candidate for the HH 7-11 outflow (Rodriguez, Anglada, &: Curiel 1997). However, source A2 is 
only a 3a detection in panel (c) and is not detected at higher resolution. We suggest that its 
lack of compact structure makes it a less likely candidate for driving the outflow, despite the 
fact that the centimeter emission suggests that it also has a jet. In addition, the high spectral 
index in the centimeter (a ~ 1.5; Rodriguez, Anglada, & Curiel 1997) and the lack of detection in 
high resolution A = 1.3 mm observations (Bachiller et al. 1998) suggest that this source may be 
dominated by free- free emission. Source Al is coincident with the infrared/optical source SVS13. 
Since source Al is an optical source, we would expect it to be an older object. However, our data 
suggest that Al is deeply embedded in an envelope. The SVS13 results are discussed in detail in 
Welch, Looney, & Mundy (1999). 



3.10. NGC 1333 IRAS4 

One of the most well known sources in the NGC 1333 region is the object NGC 1333 IRAS4, 
located to the southwest of SVS13. Unresolved in the IRAS images (Jennings et al. 1987), 
IRAS4 breaks into two bright objects at sub-millimeter wavelengths (Sandell et al. 1991). Our 
images. Figures 12, 13, and 14, show three objects: IRAS4 A, B, and C. Our data provide the 
first indication that source C may be a young star. Source C is detected at all resolutions, has a 
5kA/50kA ratio near 1, and the integrated flux in Table 2 is constant at all resolutions. The source 
C characteristics are more like those of an optical/IR source than its IRAS4 companions. 

NGC1333 IRAS4 A & B have been observed with the CSO-JCMT single basehne 
interferometer at A = 840 fini (Lay, Carlstrom, 8z Hills 1995). Their best fit for source A was two 
elliptical Gaussians, and indeed, in our images source A is a binary system. It is interesting to 
note that the best fit fluxes from Lay et al. give a ratio of 0.78, while our two sources have a flux 
ratio of 0.25. This suggests that either the dust emissivity of these two objects vary differently 
with frequency or the optical depth is very different. For source B, the CSO-JCMT data could not 
be fit with a single star or binary model. Lay et al. suggest that source B may be a triple system; 
however, they were not aware of source C at that time, which may have confused their analysis. 
Our image of source B shows weak extensions (4(j) to the north and southwest, but our data are 
not sufficient to determine the nature of these features. They could trace a multiple stellar system 
or inhomogcncitics within the envelope. 
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3.11. VLA 1623 

The source VLA 1623, near the center of the p Ophiuchi cloud core A, is the prototypical 
Class source (Andre, Ward-Thompson, &; Barsony 1993) that drives a large outflow with 
principal axis ~ -60° (Andre et al. 1990; Dent et al. 1995; Yu & Chernin 1997). This source 
has been observed with the CSO-JCMT single baseline interferometer at A = 1360 and 845 /xm 
(Pudritz et al. 1996). They modeled the source as a Gaussian and placed a 70 AU radius upper 
limit on the size of the compact circumstellar disk. 

Recent, high resolution VLA observations at A = 3.6 cm (Bontemps & Andre 1997) show a 
series of emission clumps that are interpreted as knots of a radio jet driving the large CO outflow. 
However, the position angle of the radio jet and the CO outflow differ by ~ 30°. In our highest 
resolution images. Figure 15 panels (c) and (d), the millimeter emission breaks into nearly equal 
point sources. The two crosses mark the positions of the two centimeter sources from Bontemps & 
Andre (1997) that appear associated with the millimeter emission; the source positions at the two 
wavelengths agree to within the uncertainties. The total emission from the two sources at A = 3.6 
cm is less than 1 mJy, so free-free emission is probably not contributing much of the 70 mJy seen 
at A = 2.7 mm. In addition, we have reanalyzed the data of Pudritz et al. (1996) and find that 
a binary interpretation is also consistent with their data due to their limited u,v sampling. Thus, 
the emission at A = 2.7 mm is probably dominated by dust emission. VLA 1623 is most likely a 
very young binary system with two circumstellar disks. Like IRAS 16293-2422, we refer to the 
southern source as A and the northern source as B. 



3.12. IRAS 16293-2422 

IRAS 16293-2422 is a very well studied deeply embedded binary system with two molecular 
outflows (Walker et al. 1986; Wootten 1989; Mundy et al. 1992; Walker et al. 1993) in p Ophiuchi. 
The outflow from the southern source A has a principal axis of ~ 50°, and the outflow of the 
northern source B has a principal axis of ~ 75°. The outflow from source B does not extend down 
near to the source, which may indicate that source B is no longer driving its outflow. In high 
resolution observations at A = 2 cm, the system is comprised of three peaks: Al and A2 to the 
southeast and B to the northwest (Wootten 1989). In Figure 16, we detect the two sources, A and 
B, that were detected previously at A = 2.7 mm (Mundy ct al. 1986; Mundy ct al. 1992). In panel 
(c), there is still a clear connection between the two sources that is most likely a circumbinary 
envelope. In panel (d), the massive circumbinary envelope is mostly resolved-out and the residual 
emission arises from two compact sources and some weak extensions that are probably attributed 
to density enhancements within the circumbinary structure. At high resolution, source A appears 
elongated along the position angle of the A = 2 cm sources, which arc indicated in panel (d) as 
crosses. IRAS 16293-2422 source A has the highest 5kA/50kA ratio in the survey. In fact, the ratio 
is twice as large as the next highest source L1448 IRS3 B. Thus, source A has most of its mass in 
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the envelope, perhaps making it one of the younger sources in this survey. 

Our measurement of the integrated flux from this source is higher than previous observations 
(e.g. Walker et al. 1993). This is because we have shorter spacing u,v data which pick up the 
extended structure of the circumbinary envelope better than previous works. If we remove the 
shorter u,v spacings, the total integrated flux is ~750 mJy, which is more in agreement with the 
previous measurements. 

4. Comparison of Structure 

There is a striking difference between the embedded objects and the optical/IR objects in our 
survey. The optical sources have compact emission on spatial scales of ~1" with little large-scale 
envelope emission. This is illustrated both in Figure 17 and by the peak/integrated fluxes in 
Table 2. The peak flux does not change significantly, even down to size scales of ~1'.'5, until the 
resolution is sufficient to see the circumstellar disk. This contrasts strongly with the embedded 
objects which typically have > 85% of their emission in large scale structures. The embedded 
sources, Figures 7 through 16, show dramatic change in structure as the resolution is varied 
through the panels. Structures are resolved-out as the shorter u,v spacings are down- weighted in 
the higher resolution images. At the highest resolution, the embedded objects typically have a 
residual compact component, but the flux of this component is significantly less than that in large 
scale extended emission. In addition, the embedded objects often show complex sub-structure 
within the field. 

Indeed, it is usually difficult to isolate the circumstellar disk from the envelope in the 
embedded sources, even at (f.'b resolution. In most cases, the emission shows no discontinuity in 
flux between 1" and 3" scales; this indicates that any disk present cannot be significantly more 
massive than the mass of the envelope extended to small scales. Since the embedded sources are 
typically a factor of two farther away than the optical sources, we cannot establish whether the 
disks in embedded systems are systematically less massive than typical disks around young optical 
stars. Nonetheless, the younger disks do not appear to be prominent mass reservoirs compared to 
their envelopes on the few hundred AU scale. 

However, it is interesting to note that even in the two nearest embedded systems VLA 1623 
and IRAS 16293-2422, where our spatial resolution is comparable to the Taurus systems, the disk 
and envelope contributions to the continuum emission are difficult to disentangle. In the VLA 
1623 system, roughly one-third of the total flux is in the large-scale circumbinary envelope, and 
two-thirds is in the two unresolved components (< 160 diameter); each compact component may 
be a ~O.O4M0 circumstellar disk or a combination of an individual circumstellar envelope and 
a circumstellar disk. In the IRAS 16293-2422 system, most of the flux (~85%) is in the large 
scale circumbinary envelope and any disk flux is less than 15% of the total flux, but even here the 
envelope emission could still be significant fraction of the flux at small scales. 



-14- 



Overall, our results generally support a picture in which most embedded systems have disk 
masses which are comparable to, or less than, the typical disk masses of ~0.02 Mq foTind for 
young optical stars (Osterloh & Beckwith 1995). Most theoretical works find that the disk grows in 
prominence with time as the system evolves. Cassen &; Moosman (1981) showed that the detailed 
evolution of a disk is very dependent upon the distribution of mass and angular momentum in the 
original cloud and dissipative processes within the disk. For a range of assumptions, they found 
that the disk should grow more massive and larger with time. Building upon these results, Stahler 
et al. (1994) considered a disk with negligible viscosity that was formed as soon as the angular 
momentum in the infalling material caused it to "miss" the protostar. They found that the radius 
of the disk is a strong function of time, increasing as t^. As the accretion rate onto the star begins 
to fall off, the mass of the disk increases. Recent results using magnetic collapse models (Basu 
1998), predict that the early disk grows less strongly, increasing only linearly with time, but, in 
these models, the disk starts out more massive since less material can fall directly onto the star. 
Unfortunately, even higher linear resolution observations than those presented here are needed to 
fully isolate disks in most embedded systems and test detailed theoretical predictions. 



5. Simple Mass Comparison 

How does the circumstellar mass in the optical systems and embedded systems compare? 
The A = 2.7 mm emission provides a measure of circumstellar mass in these systems. Detailed 
modeling, including data from a variety of wavelengths, is needed to provide the best mass 
estimates, but the accuracy of any mass estimate is limited by the uncertainty in the dust 
emissivity and the true structure of the source. Since the purpose of this paper is to present a 
survey, we will apply a very simple emissivity model, and some assumptions, to make a rough mass 
estimate for the different systems. We adopt a single-temperature, optically thin source model 
so that Fjy = B^{Tiiust)i<-uM/D^ , where B^{T) is the Planck function, T^ust is the temperature of 
the dust, K,y is the dust mass opacity, M is the mass of gas and dust, and D is the distance to 
the source. For the dust temperature, we will assume characteristic temperatures of 35 K for the 
deeply embedded objects and 60 K for the disk dominated systems. For the dust mass opacity 
(Beckwith &; Sargent 1991; Draine 1990; Pollack ct al. 1994; Stognienko, Henning, & Ossenkopf 
1995), we adopt a which is consistent with other works (e.g. Beckwith & Sargent 1991; Ohashi 
et al. 1991; Osterloh & Beckwith 1995): = 0.1(i//1200 GHz) cm^ g^^, corresponding to Ki, = 
0.009 cm^ g-i at A = 2.7 mm. Although we do not expect this simple model to give accurate 
masses, it provides rough estimates that are adequate for qualitative comparisons and are arguably 
within a factor of 2 of the likely mass. More detailed modeling of the individual sources will be 
done in subsequent papers. 

Table 3 lists the estimated mass for each source, as well as the best fitted positions from the 
highest resolution image (typical positional uncertainties of 0'.'15). The masses show more than a 
factor of a hundred range from the least to most massive system in our sample. Where there is 
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overlap, there is reasonable agreement between the simple model and published mass estimates. 
For example, the mass for HL Tauri in Table 3, 0.05 Mq, is within the range of masses previously 
found, 0.05 to 0.1 Mq (Beckwith et al. 1990; Mundy et al. 1996, Wilner, Ho, & Rodriguez 1996; 
Close et al. 1997), and our mass for DG Tauri, 0.03 Mq, is consistent with previous estimates of 
0.02 to 0.04 M0 (Beckwith et al. 1990; Dutrey et al. 1996). For the embedded systems IRAS 
16293-2422 and L1448 IRS3, our estimated masses of 1.1 Mq and 0.74 Mq respectively, are 
similar to other interferometric estimates, 1 Mq for IRAS 16293-2422 (e.g. Mundy et al. 1992) 
and 0.5 Mq for L1448 IRS3 B (e.g. Terebey, Chandler, & Andre 1993). 

The circumstellar masses for the two categories, optical/infrared and deeply embedded 
sources, follow the expected trend: the embedded objects typically have a factor of 5 or so larger 
masses. Comparing circumstellar masses and stellar masses for the optical sources, as expected, 
these stars have already acquired most of their final mass. The luminosities of the optical/infrared 
systems range from ~1 to 30 Lq, suggesting central masses of 0.5 to 1.5 Mq, whereas their 
circumstellar masses in Table 3 range from 0.01 to 0.08 Mq. 

The embedded systems have typical circumstellar masses of ~0.5 Mq, with the largest one, 
NGC 1333 IRAS 4 A, near 3 Mq. The luminosities of the embedded systems range from ~ 1 Lq 
for VLA 1623 to ~ 50 Lq for the NGC 1333 SVS 13 system. Given the star formation regions 
in which they are found (NGC 1333 and Ophiuchi) and their luminosities, it is likely that the 
embedded sources are forming a similar range of stellar masses to that of the optical/infrared 
sources. Thus, the current circumstellar masses are generally comparable to the expected stellar 
masses. Unfortunately, the current masses of the forming stellar bodies is unknown. Nearly all 
of the observed systems have outflows, centralized dust heating, and significant luminosities, 
indicating that a stellar object is present. The current stellar mass required to generate the 
observed luminosity can be roughly estimated, if it is assumed that the observed circumstellar 
mass is accreting on a central object at a constant rate over the typical duration of this phase of 
10^ years. For a stellar radius of 3 Rq, the derived stellar masses are in the range of 0.05 Mq to 
0.4 Mq. This simple approach suggests that the stellar cores in the embedded systems are likely 
to have already attained a significant fraction (10 to 40%) of their final stellar masses. 

6. Young Multiple Systems 

All of the embedded sources in our survey are either part of small groupings or are in close 
binary systems. Even though our sample may be biased due to our selection criteria, binary 

systems appear to be common in the earliest stages of star formation. Three possible explanations 
for the large number of embedded binary systems are: (1) small sample statistics, (2) selection 
effects in our sample, such that binary, embedded sources were preferentially chosen, or (3) true 
high multiplicity among younger systems (Ghez, Neugebaucr, &; Matthews 1993). 

On the first point, even though we observed a small number of systems, the probability of 
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choosing six multiple systems from a Taurus-like binary fraction distribution of ~ 60% (e.g. Ghez, 
Neugebauer, &, Matthews 1993) is low (only 4.7%) but it can not be ruled out. On the second 
point, one obvious selection bias would be if binary systems typically have more massive envelopes 
than single star systems, making embedded binary systems brighter at millimeter wavelengths. 
This supposition is opposite to the trend seen in older, optical T Tauri binaries; studies of these 
systems (Beckwith et al. 1990; Jensen, Mathieu, & Puller 1994, 1996; Osterloh & Beckwith 1995) 
provide statistical evidence that T Tauri binary systems have less millimeter emission than single 
systems. However, the emission properties of young binary systems could be time dependent — 
young embedded binary systems could be brighter millimeter sources than coeval single star 
systems, because they have more massive envelopes. As they evolve, their envelopes disappear 
and the remaining material is more rapidly processed through their circumstellar disks than in 
single star systems, and they become less bright at millimeter wavelengths than comparable single 
star systems. Data on more embedded systems are needed to test this possibility. 

On the third point, the evolution of binarity among stellar systems is an open question. Our 
results support the idea that binaries and multiple systems may be more common in younger 

systems; however, our criteria for identifying multiple systems is fairly loose. Systems with 
separations of 5000 AU or more arc valid pairs for forming stars because the reservoir of cloud 
material that the forming stars draws from is typically several thousand, or more, AU. But, as 
such systems evolve, the loss of envelope mass and interactions with other stars forming in the 
cloud provide mechanisms for unbinding loose binary systems. The evidence for evolutionary 
trends in the binarity of optical young stars is contradictory. Speckle observations of the Hyades 
cluster, a young main-sequence cluster, find that the occurrence of binary systems in the cluster 
is larger than in the local solar neighborhood but less than the Taurus clouds (Patience et al. 
1998), but studies of additional optical clusters did not confirm the trend of decreasing binary 
occurrence with age (Patience et al. 1999). It is also possible that binary formation mechanisms 
have a dependence on the initial conditions of the pre-collapse cloud (Durisen Sz Sterzik 1994) or 
the density of the star forming cluster (Bouvier, Rigaut, &; Nadeau 1997), so simple comparison 
of different regions may not be valid. Broader survey work is needed to separate the different 
possible factors. The solid result from our work is that multiplicity is common and established 
early in the formation process and on a variety of scales. 

The most favored mechanism for the early formation of binary and multiple stellar systems is 
fragmentation within either the initial cloud core or the circumstellar disk. Fragmentation during 
the earliest stages of the isothermal collapse of a cloud core, due to perturbations or non-spherical 
cores, can form binary systems with separations ranging from 10 to 10^ AU (Boss &; Bodenheimer 
1979; Monaghan & Lattanzio 1986; Bonnell et al. 1991; Bonncll & Bastien 1992; Boss 1993; 
Bate, Bonnell, & Price 1995). Fragmentation due to gravitational instabilities in the circumstellar 
disk systems may form binary systems with separations ranging from 10 Rq to 100 AU (Adams, 
Ruden, & Shu 1989; Shu et al. 1990; Bonneh 1994; Bonnell & Bate 1994). In our survey, the 
majority of the circumstellar mass in embedded systems is in the large-scale envelope, with very 
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little mass in circumstellar disks. This would suggest that fragmentation occurs during the early 
evolution of the core in most systems. However, the Perseus clouds are distant enough that we 
would not detect close (< 150 AU separation) binary systems, the primary regime of the disk 
fragmentation and the dominate population for optical binary stars. Of the optical sources in our 
survey, all of which are located in the Taurus cloud, only two appear to be binary systems — GG 
Tauri and L1551 IRS5. Both of these sources could have been formed by disk fragmentation since 
they both have circumbinary structures. Observations with resolution of 10 AU, or smaller, are 
needed to probe the fragmentational history of close binaries. 

Morphologically, we can identify three types of multiple systems in our sample: separate 
envelope, common envelope, and common disk systems (Table 4). The characteristics of the 
different systems are defined by the broad distribution of the circumstellar material. Separate 
envelope systems exhibit clearly distinct centers of gravitational concentration with separations of 
>6000 AU; the components are within a larger surrounding core of low density material. Common 
envelope systems have one primary core of gravitational concentration which breaks into multiple 
objects at separations of 100 - 3000 AU. Common disk systems have separations of < 100 AU and 
typically have circumbinary disk-like distributions of material. Table 5 lists the binary systems 
with our classification, their association, and the projected separation. The association number 
assigned in Table 5 identifies members of common envelope or common disk systems. 

There are several clear connections between these morphological distinctions and other works. 
The study of the separation distribution of optical binaries by Larson (1995) found a knee in the 
distribution at 0.04 pc (8250 AU) which was identified with the Jeans size. Larson suggested 
that systems on that scale and larger formed by fragmentation and separate collapse, exactly the 
structure found in the separate envelope systems. This scenario of prompt initial fragmentation 
is not new (e.g. Larson 1978, Pringle 1989, Bonnell et al. 1991); it was discussed recently by 
Bonnell et al. (1997) in the context of small cluster formation. The critical issue is that the 
collapse is initiated in a system which contains multiple Jeans masses in a weakly condensed 
configuration; one example of such a system would be a prolate Gaussian distribution with several 
Jeans masses along the long axis and one Jeans mass across the short axes. In systems with large 
separations, ambipolar diffusion in self-gravitating, magnetic clouds may play an important role 
(e.g. Mouschovias 1991). 

The common envelope systems can be linked with models for the fragmentation of moderately 
centrally-condensed spherical systems (Burkert & Bodenheimer 1993; Boss 1995, Boss 1997). In 
these models, fragmentation occurs in the dense central region within an overall single core. The 
primary requirement for fragmentation is that the central region has a fairly flat distribution. 
However, evidence of this flat region is erased quickly once the fragmentation and collapse occurs, 
and the forming multiple system is left embedded within a single centrally condensed core. Finally, 
the common disk systems are similar to models of high angular momenta systems (Artymowiez 
& Lubow 1994; Bate & Bonnell 1997). This implies that close stellar systems may be created 
by fragmentation of early disks or pseudo-disk structures created early in the collapse (e.g. Galli 
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& Shu 1993). The distribution of material between circumstellar and circumbinary structures in 
such systems, depends sensitively on the angular momentum of the infalling material. 

7. Conclusions 

We have presented the first sub-arcsecond millimeter wavelength survey of the dust continuum 
emission towards 24 young stellar systems. Morphologically, the systems can be divided into 
young, embedded objects and older optical/infrared sources. The optical systems show compact 
emission from circumstellar disks that only begins to be resolved at ~1" resolution. In two systems, 
HL Tauri and DG Tauri, we resolve the circumstellar disk. The embedded systems show thermal 
dust emission that is dominated by their circumstellar envelopes with weak residual emission at 
small scales. The embedded systems have > 85% of their A = 2.7 mm continuum emission from 
large scale structures. This suggests that young circumstellar disks of the embedded systems are 
not prominent mass reservoirs compared to the inner circumstellar envelope. In fact, for most of 
the embedded systems, there is no obvious discontinuity in flux between 1" and 3". Thus, if there 
are circumstellar disks in these systems, they are not more massive than the expected mass of the 
power-law envelope extrapolated to small scales. 

Simple estimates of the circumstellar mass (not including the star) in our systems show that 
the circumstellar mass of the embedded systems (ranging from 0.04 to 2.9 Mq) is, as expected, 

larger than the circumstellar mass of the optical/infrared systems (0.01 to 0.08 Mq). Typically, 
the embedded systems have a factor of 5 or so larger circumstellar masses. This follows the broad 
trend that the older optical/infrared sources have nearly reached their final stellar mass, since 
there is little mass left to accrete, and that the younger, embedded sources have not reached their 
final stellar mass, since it is expected that a fraction of the circumstellar mass will be accreted. 
However, simple calculations for the embedded systems based on luminosity and circumstellar 
mass, suggest that the embedded systems may have already accreted a significant fraction (10% 
or more) of their probable final stellar mass. 

The survey has a large number of multiple systems; all of the embedded objects are in small 
groupings or binary systems. Morphologically, based on the system separation, we categorize our 
sample into three types of multiple systems: separate envelope (separation > 6500 AU), common 
envelope (separation 150-3000 AU), and common disk (separation < 100 AU). The separate 
envelope group links with the idea of initial fragmentation and collapse of multiple Jeans masses 
within a weakly condensed configuration (e.g. Larson 1978; Pringle 1989; Bonnell et al. 1991). The 
common envelope group is similar to fragmentation of a moderately centrally condensed spherical 
system (e.g. Burkcrt & Bodenheimer 1993; Boss 1995, 1997). Finally, the common disk group 
connects with fragmentation of the early disk due to high angular momentum (Artymowiez & 
Lubow 1994; Bate k, Bonnell 1997). These three categories suggest three distinct times during the 
star formation process where specific fragmentation mechanisms dominate the binary formation 
process. Thus, the separation morphology may also indicate a time sequence at which binaries 
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form in stellar evolution. 
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Table 1. Source List 



Source 


Distance 
(pc) 


Optical/IR 
or Embedded 


Class 


Main 
Calibrator 


Secondary 
Calibrator 


Dist. 
Ref. 


L1448 IRS3 A 


300 


Embedded 





3C111 


0336+323 


2 


L1448 IRS3 B 


300 


Embedded 





3C111 


0336+323 


2 


L1448 IRS3 C 


300 


Embedded 





3C111 


0336+323 


2 


NGC1333 IRAS2 A 


350 


Embedded 





3C111 


0336+323 


2 


NGC1333 IRAS2 B 


350 


Embedded 





3C111 


0336+323 


2 


SVS 13 Al 


350 


Optical/IR 




3C111 


0336+323 


2 


SVS 13 A2 


350 


Embedded 




3C111 


0336+323 


2 


SVS 13 B 


350 


Embedded 




3C111 


0336+323 


2 


SVS 13 C 


350 


Embedded 




3C111 


0336+323 


2 


NGC1333 IRAS4 Al 


350 


Embedded 





3C111 


0336+323 


2 


NGC1333 IRAS4 A2 


350 


Embedded 





3C111 


0336+323 


2 


NGC1333 IRAS4 B 


350 


Embedded 





3C111 


0336+323 


2 


NGC1333 IRAS4 C 


350 


Embedded 




3C111 


0336+323 


2 


DG Tauri 


140 


Optical/IR 


II 


3C111 


0431+206 


1 


DG Tauri B ^ 


140 


Optical/IR 


I 


3C111 


0431+206 


1 


L1551 IRS5 


140 


Optical/IR 


I 


3C111 


0336+323 


1 


HL Tauri 


140 


Optical/IR 


II 


0530+135 


0431+206 


1 


GG Tauri 


140 


Optical/IR 


II 


0530+135 


0431+206 


1 


GM Aurigae 


140 


Optical/IR 


II 


3C111 


3C123 


1 


VLA 1623 


160 


Embedded 





1733-130 


1625-254 


3 


IRAS 16293-2422 A 


160 


Embedded 





1733-130 


1625-254 


3 


IRAS 16293-2422 B 


160 


Embedded 





1733-130 


1625-254 


3 



^DG Tauri B observed near the FWHM of primary beam. 

References. — (1) Elias 1978; (2) Herbig k Jones 1983; (3) Whittet 1974 
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Table 2. Source Flux 



OUUI 




JrcdK S: ILIA 

y Vl 1 T r / r\£iOTVl 

ilijy / Ucdiiii 


iiibcgiaitjo p iux 
mJy 


DUX. Oizlc 


L1448 IRS3 A 


(a) 


26.5±1.6 


23.1± 2.6 


11^'2 X 6^'3 




(b) 


14.5±1.5 


26.7± 3.3 


8^'3 X 5'.'8 




(c) 

\ / 


6.8±1.6 


19.3± 3.7 


2^'3 X 2^'9 




(d) 


<6.7 






L1448 IRS3 B 


(a) 


101.5±1.6 


134.6± 3.9 


17'.'0 X 9'.'8 




(b) 


84.9±1.5 


135.6± 4.8 


ll'.'O X g'.'o 




(c) 

\ / 


41.1±1.6 


135.2± 6.7 


4?5 X 4'.'9 




(d) 


22.5±2.3 


115.7± 9.5 


2:'8 X 2'.'4 


L1448 IRS3 C 


(a) 


14.9±1.6 


31. 7± 4.1 


lQf.'7 X 17'.'0 




(b) 


11.7±1.5 


31.9± 3.7 


S'.'O X 12'.'0 




(c) 


9.8±1.6 


14.3± 3.2 


2'.'5 X 2'.'0 




(d) 


8.7±2.3 


8.7± 2.3 


0^'7 X 0'.'5 


N1333 IRAS2A 


(a) 


46.5±1.3 


82.8± 4.0 


16'.'4 X 16?3 




(b) 


36.2±1.2 


74.4± 4.0 


ll'.'O X 12'.'0 




(c) 


22.3±1.7 


36.1± 4.4 


2'.'9 X 2'.'2 




(d) 


18.4±2.7 


22.4± 4.8 


1'.'4 X o'.'g 


N1333 IRAS2B 


(a) 


21.3±1.3 


27.7± 3.2 


12'.'8 X 13'.'0 




(b) 


19.6±1.2 


24.4± 2.7 


6?6 X 9'.'2 




(c) 


18.9±1.7 


24.7± 3.5 


1^'9 X 2'.'2 




(d) 


16.9±2.7 


24.3± 5.1 


1^'2 X 1'.'2 


SVS 13 A 


(a) 


47.3±1.1 


101. 3± 4.2 


26'.'0 X IS'.'O 




(b) 


37.4±1.1 


100.3± 4.7 


20'.'7 X 10'.'2 




(c) 


19.2±1.5 


38.7± 4.1 


2'.'9 X 2'.'9 




(d) 


11.0±2.2 


38.0± 6.6 


1'.'7 X 2'.'2 


SVS 13 B 


(a) 


52.0±1.1 


123.0± 4.5 


31'.'5 X 14':'0 




(b) 


41.7±1.1 


110.4± 3.7 


10'.' 1 X 13'.'0 




(c) 


25.3±1.5 


41.4± 3.6 


2'.'2 X 3'.'0 




(d) 


19.4±2.2 


48.2± 6.6 


I'.'S X 2'.'5 


SVS 13 C 


(a) 


11.7±1.1 


21. 0± 2.5 


14'.'0 X lO'.'O 




(b) 


9.6±1.1 


19.8± 2.7 


9^'l X 7'.'6 




(c) 


8.5±1.5 


8.5± 1.5 


I'.'l X I'.'O 




(d) 


11.1±2.2 


11. 1± 2.2 


0'.'7 X O'.'S 


N1333 IRAS4A 


(a) 


351.2±3.1 


544.2±13.6 


25'.'0 X 24'.'0 




(b) 


280.4±1.9 


525. 6± 9.2 


12'.'0 X 18'.'5 




(c) 


172.2±2.1 


449.7± 9.8 


5'.'4 X 6'.'2 


Al Only 


(d) 


107.0±2.9 


324.1±12.0 


2^'9 X 2'.'2 
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Table 2 — Continued 



Source Panel Peak Flux Integrated Flux Box Size 
mJy/beam mJy 



A2 Only (d) 



23.0±2.9 



81.2± 8.1 



I'.'S X l'/6 



-28- 



Table 2 — Continued 



Source 


Panel 


Peak Flux 

T /I 

mJy/beam 


Integrated Flux 
mJy 


Box Size 


N1 '^'^'^ TR A^4R 


(a) 




1 on 7 q 


1 9"0 V 1 7"n 




(h) 


1 oq 1 +1 q 


1 72 1 + n 


8"5 X 1 1"0 




(c) 


94.0±2.1 


148 Q± 5 Q 


3'.'4 X 3''6 








1 90 o_|_ 7 q 


V'7 X T'fi 


N1 TR AS4C 


fa) 


47 Q + Q 1 


4q 01 5 5 


q"n X 1 i"o 






"-±0. <J_I — L . U 


51-) 7_)_ Q 


5"5 X 7"0 




(c) 


Qq q+2 1 


57 n_|_ 4 q 

'~J I .\J—\— '-t.U 


3"0 X 2"8 






na q_Lo q 


B1 6+ 6 Q 


1"4 X 1"5 

X . ^ /\ X . (J 


DG Tauri 

J—' VJI _L Ct LI J- -L 




57.7±2.7 


66 0± 5 8 


11''4 X 11'.'2 






fi+2 f) 


70 01 e 9 


Q"0 X 10"4 




(c) 




71 9+ 4 


2"q X 2"9 








Ro qi 5 1 


9"l X 9"5 

. X /S ^ , ' J 


DO Taiiri R 




45 n+4 Q 


70 4_|_-| 1 Q 


1 3"0 X 1 2"0 








79 7-1-1 n q 


7"8 X 11 "8 

1.0 r\ X X . 




(c) 




47 o-L 7 9 

t: f . 0_1_ 1 


1'.'9 X 3''1 




I'd) 




4q 4J- Q 


VIA y -i/zo 

X . ^ /\ X . 




[A) 




1 70 O-L 7 5 


17"0 X 14"0 




(h) 


1 on 7+0 5 


1 77 O-L 7 q 


12"0 X 8"8 

XZ. VJ /V (J, 




(c) 


77 q_|_Q Q 


145 O-L q 1 


9" 7 X "^''^ 

I /N (J. U 




(d) 




107 flill 1 

X f • 1 L J. • J. 


1'.'3 X 1''9 


HT, Tfliiri 


yd.) 


102 7+1 7 

X.\j£j. % _l L . 1 


1 08 6+ 4 6 


1Q"0 X 17"8 

Xc/.L/ /N XI. 






qn q_|_-| 7 


1 1 fJ-L 4 


Q"0 X Q"fi 






7n Q_|_9 4 


-|n« 9-L « n 


9"fi X 9"q 




Cd) 


4a c+9 q 


-|n« q+ 7 


I'.'S X I'.'S 


CiG Tauri 


yd) 


5« 7+1 

(JLI. 1 —1 — L .0 


7Q 5-L 4 4 


1 2"6 X 1 2"4 




(h) 


27.2±1.2 


72 5± 3 6 


7% X 6''7 




(c) 


10.3=bl.2 


78 0± 5 


4'.'8 X 4''9 




(d) 


10.8±1.5 


95.2± 6.8 


4'.'4 X 4'.'5 


GM Aurigae 


(a) 


20.3±1.1 


22.0± 2.6 


12'.'0 X 14'.'0 




(b) 


19.2±0.9 


22. 3± 2.0 


7^8 X 6'.'2 




(c) 


13.6±1.6 


19.6± 3.0 


2'.'0 X 2'.'2 




(d) 


13.4±2.5 


13.4± 2.5 


0'.'6 X 0^'5 


VLA 1623 A&B 


(a) 


54.2±3.0 


72.1± 6.8 


9^'4 X 18'.'0 




(b) 


44.2±2.2 


53.5± 4.3 


5'.'4 X 8'.'0 


A Only 


(c) 


22.8±2.0 


34.4± 4.3 


1'.'6 X 3'.'4 


B Only 


(c) 


25.0±2.0 


32.5± 4.3 


1'.'6 X 3'.'4 
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Table 2 — Continued 



Source Panel Peak Flux Integrated Flux Box Size 
mJy/beam mJy 

A Only (d) 22.4±3.5 25.5± 6.3 (y/9 x l'/7 

B Only (d) 25.8±3.5 25.8± 3.5 (X'4 x Qf.'9 
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Table 2 — Continued 



Source 


Panel 


Peak Flux 


Integrated Flux 


Box Size 






mJy/beam 


mJy 




IRAS 16293-2422 


(a) 


412.6±5.8 


1017.9±26.5 


22'/2 X 27^'2 


A Only 


(b) 


176.4±4.2 


441.2±14.1 


12'/! X icy/o 


B Only 


(b) 


305.3±4.2 


551.4±14.1 


12'.'5 X 9'.'7 


A Only 


(c) 


60.1±4.1 


358.3±18.5 


5'.'2 X 5'.'1 


B Only 


(c) 


154.1±4.1 


498.4±17.2 


e'.'O X 3'.'8 


A Only 


(d) 


43.6±4.8 


276.2±22.9 


3'/l X 4^8 


B Only 


(d) 


112.7±4.8 


424.2±24.2 


3'/4 X 4'/9 



'^DG Tauri B fluxes were corrected for primary beam attenuation; thus 
fluxes given have a larger overall uncertainty than the rest of the survey. 
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Table 3. Positions, Simple Estimates of Mass, and 5kA/50kA Ratio 



Source 


a (J2000) 


5 (J2000) 


Mass 


5kA/50kA 








(Mq) 


Ratio 


LI 448 TRS8 A 


n3'*25™36'^532 


+30°45'2l"35 


0.09 


3.7±0.4 


L1448 IRS3 B 


03'*25'"36'^339 


+30°45'l4"94 

1 KJ\J ^.\J -Lit t_/ A 


0.52 




L1448 IRS3 C 


03'*25™35'^653 


+30°45'34"20 


0.12 




NGC1333 IRAS2 A 


03'*28"*55^571 


-l-31°14'37"22 


0.44 


2.5±0.4 


NGC1333 IRAS2 B 


03'^28"*57^349 


-l-31°14'l5"93 


0.14 


l.libO.2 


SVS 13 Al 


03'^29"*03^750 


-l-31°16'03"95 


0.54 ^ 


3.3±0.6 ^ 


SVS 13 A2 


03'*29"*03^374 


-l-31°16'oi"87 


0.54 ^ 


3.3±0.6 ^ 


SVS 13 B 


OS'* 29"* 03^ 056 


+31°15'51"67 


0.65 


2.6±0.4 


SVS 13 C 


03'*29"*0FQ51 


-l-31°15'38"27 

1 'J J- -L (JO. ^1 


0.11 




NGC1333 IRAS4 Al 


03'*29"*lfl^510 


-l-31°13'31"fll 


2.88 ^ 


2.8±0.1 ^ 


NnC1333 TRAS4 A2 


03^*29"* 10^41 3 


Vs2"20 


2 88 ^ 


2.8±0.1 ^ 


NnC13S3 TRAS4 B 


r)o/ioqm-| isqes 


-1-31 1 3 08"! n 


96 


1 8-1-0 1 


NC;C1833 IRAS4 C 


03^29™ 12*^81 3 


-l-31°13'nfi"97 


0.26 


l.OzbO.l 


DG Taiiri 


04^27™ (14*^686 

V^t: { \j\-t\J\J\J 


+26° Ob' 16" 14 


0.03 


0.9±0.1 


UG lauri B 


U4 27 Uz.5dz 


-|-ZD U5 3U. 53 


0.04 


2.2±0.4 


L1551 IRS5 A 


04'*31™34!143 


-M8°08'05'/09 


0.08 ^ 


1.8±0.1 * 


L1551 IRS5 B 


04'' 31™ 34! 141 


-M8°08'04'.'74 


0.08 ^ 


1.8±0.1 ^ 


HL Tauri 


04'*31"^38!413 


-M8°13'57'.'61 


0.05 


1.3±0.1 


GG Tauri 


04^*32"* 30?322 


+17°3l'4a'65 


0.03 


4.0±0.6 


GM Aurigae 


04^*55"* 10?983 


+30°2l'59^'37 


0.01 


1.1±0.2 


VLA 1623 A 


16'*26"*26?396 


-24°24'30^'45 


0.04 


2.4±0.3 ^ 


VLA 1623 B 


16'*26"*26?318 


-24°24'30'/12 


0.04 


2.4±0.3 ^ 


IRAS 16293-2422 A 


16'*32"*22!869 


-24° 28 36'/ 11 


0.49 


7.4±0.4 


IRAS 16293-2422 B 


16'*32"*22!624 


-24°28 32'.'20 


0.61 


3.5±0.1 



^Closc binary systems whose values include both systems. 

^DG Tauri B observed at FWHM point of primary beam; thus masses 
given have a larger uncertainty than the rest of the survey. 

Note: The mass scales roughly linear with the assumed temperature. For 
example, in the case of the embedded objects, M(T) = M(35K)x(35K/T). 
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Table 4. Characteristic Scales for Multiplicity 



Classification 


Scale (AU) 


Separate Envelope 


> 6000 


Common Envelope 


150 - 3000 


Common Disk 


< 100 
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Table 5. Multiple System Morphology 



Source 


Type 


Assoc. 


Separation 










Arcsec 


AU 


L1448 IRS3 A 


common 


envelope 


1 


D. 87 


z4U4 


L1448 IRS3 B 


common 


envelope 


1 


aft 0*7 
D. 87 


z4U4 


L1448 IRS3 C 


separate 


envelope 




1 '7// 1 O 

17. 13 


5995 


NGC1333 IRAS2 A 


separate 


envelope 




Ol. zU 


1 nnon 
lU9zU 


NGC1333 IRAS2 B 


separate 


envelope 




ol. 10 


1 nnon 
lU9zU 


SVS 13 Al 


common 


envelope 


2,3 


5'.'25 


1838 


SVS 13 A2 


common 


envelope 


2,3 


5'.'25 


1838 


SVS 13 B 


common 


envelope 


3 


1(X'98 


3843 


SVS 13 C 


separate 


envelope 




19'.'50 


6825 


NGC1333 IRAS4 Al 


common 


envelope 


4 


1'.'72 


602 


NGC1333 IRAS4 A2 


common 


envelope 


4 


1'.'72 


602 


NGC1333 IRAS4 B 


separate 


envelope 




29^74 


10409 


NGC1333 IRAS4 C 


separate 


envelope 




1(X'64 


3724 


DG Tauri 


separate 


envelope 




54'.'85 


7540 


DG Tauri B 


separate 


envelope 




54'.'85 


7540 


L1551 IRS5 A 


common 


disk 


5 


0'.'35 


49 


L1551 IRS5 B 


common 


disk 


5 


0'.'35 


49 


GG Tauri 


common 


disk 


6 


(X'25 


35 


VLA 1623 A 


common 


envelope 


7 


l^'ll 


178 


VLA 1623 B 


common 


envelope 


7 


I'.'ll 


178 


IRAS 16293-2422 A 


common 


envelope 


8 


5'.' 14 


822 


IRAS 16293-2422 B 


common 


envelope 


8 


5'.' 14 


822 



^Column 3 indicates members of a multiple system. 
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DG Tau A = 2.7 mm Continuum Emission 




4''27°'4f8 4'6 4''27"°5=0 4f5 



Right Ascension (J2000) 

Fig. 1. — DG Tauri maps of the A = 2.7 mm continuum emission. All panels are contoured in steps 
of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28) x a rms noise of 2.0 mJy/bcani. (a) a = 2.7 mJy/beam; 
beam is 5'.'37 x 4'.'57 P. A. = 72°. (b) a = 2.0 mJy/beam; beam is 3'.' 12 x 2'.'72 P. A. = 68°. (c) a 
= 1.9 mJy/beam; beam is 1'.'12 x l'.'02 P.A. = 45°. (d) a = 1.6 mJy/beam; beam is 0'.'76 x 0'.'58 
P. A. = 56°. The cross in panel (d) is the A = 6 cm peak from Bieging, Cohen, & Schwartz (1984). 
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DG Tauri B A = 2.7 mm Continuum Emission 




Right Ascension (J2000) 

Fig. 2. — DG Tauri B maps of the A = 2.7 mm continuum emission. All panels are contoured 
in steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28) x a rms noise of 2.0 mJy/beam. (a) a = 2.7 
mJy/beam; beam is 5^'37 x 4('57 P.A. = 72°. (b) a = 2.0 mJy/beam; beam is 3^'12 x ^!72 P.A. 
= 68°. (c) a = 1.9 mJy/beam; beam is 1'.'12 x l'.'02 P.A. = 45°. (d) a = 1.6 mJy/beam; beam is 
0'.'76 x O'.'SS P.A. = 56°. The cross in panel (d) is the A = 3.6 cm peak from Rodriguez, Anglada, 
& Raga (1995). 
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Fig. 3. — L1551 IRS5 maps of the A = 2.7 mm continuum emission. All panels are contoured 
in steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28) x a rms noise of 3.9 mJy/beam. (a) a = 2.6 
mJy/beam; beam is 5'.'15 x 5'.'05 P.A. = -62°. (b) a = 2.5 mJy/beam; beam is 3('13 x 2'.'92 P.A. 
= 31°. (c) a = 3.3 mJy/beam; beam is l^'ll x O'.'SS P.A. = 62°. (d) a = 3.9 mJy/beam; beam is 
0'.'74 X 0'.'36 P.A. = 46°. The two crosses in panel (d) are the A = 1.3 cm peaks from Koerner & 
Sargent (1999). 
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HL Tau A = 2.7 mm Continuum Emission 



o 
o 
o 

CV2 



18 1345 



o 
B 
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4''31'°39= 




14 00 



18 13 55 




4''31°'3B?5 38=4 36=3 4''31'"38=6 38?4 

Right Ascension (J2000) 



38f2 



Fig. 4. — HL Tauri maps of the A = 2.7 mm continuum emission. All panels are contoured in steps 
of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28) x a rms noise of 2.9 mJy/beam. (a) a = 1.7 mJy/beam; 
beam is 5'.'31 x 4'.'79 P.A. = -81°. (b) a = 1.7 mJy/bcam; beam is 3'.'43 x 2'.'79 P.A. = 1°. (c) a 
= 2.4 mJy/beam; beam is I'.'ll x 0'.'94 P.A. = 53°. (d) a = 2.9 mJy/beam; beam is 0'.'68 x 0'.'48 
P.A. = 43°. The cross in panel (d) is the A = 3.6 cm peak from Rodriguez et al. (1994). 
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GG Tau X = 2.7 mm Continuum Emission 




100 AU 



4''3a'°30f4 30?2 4*32'"30?6 30?4 

Right Ascension (J2000) 



30?3 



Fig. 5. — GG Tauri maps of the A = 2.7 mm continuum emission. Panel (a) is contoured in 
steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28) x the rms of panel (a) of 1.8 mJy/beam. Panels 
(b) through (d) are (-4 -3 -2 2 3 4 5 6 8 10 14.14) x a rms noise of 1.5 mJy/beam. (a) a = 1.8 
mJy/bcam; beam is 5'.'02 x 4'.'44 P.A. = 7°. (b) a = 1.2 mJy/beam; beam is 2'.'44 x 2'.'12 P.A. = 
1°. (c) cr = 1.2 mJy/beam; beam is 1'.'17 x 1'.'02 P.A. = 31°. (d) a = 1.5 mJy/beam; beam is 1'.'02 
x 0^'80 P.A. = 40°. The greyscale is used to emphasize the hills and valleys of the "clumps" in the 
circumbinary disk. 
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GM Aurigae A = 2.7 mm Emission 
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Fig. 6. — GM Aurigae maps of the A = 2.7 mm continuum emission. All panels are contoured 

in steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28) x a rms noise of 2.5 mJy/beam. (a) a = 1.1 
mJy/beam; beam is 5'.'15 x 5'.'08 P.A. = 6°. (b) a = 0.9 mJy/bcam; beam is 3'.'12 x 3'.'00 P.A. = 
-30°. (c) a = 1.6 mJy/beam; beam is l'.'07 x I'.'OO P.A. = 51°. (d) a = 2.5 mJy/beam; beam is 
0'.'63 x 0^'47 P.A. = 62°. 
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L1448 IRS3 A = 2.7 mm Emission 




3''a5°'38?4 36f2 3''25'"38'8 36?3 



Right Ascension (J2000) 

Fig. 7. — L1448 IRS3 maps of the A = 2.7 mm continuum emission. All panels are contoured in 
steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40 56.56) x a rms noise of 2.3 mJy/beam. (a) a 
= 1.6 mJy/beam; beam is 5'.'22 x 4'.'89 P.A. = -71°. (b) a = 1.5 mJy/bcam; beam is 3'.'06 x 2'.'95 
P.A. = -61°. (c) a = 1.6 mJy/beam; beam is 1'.'08 x 0'.'99 P.A. = 56°. (d) a = 2.3 mJy/beam; 
beam is 0'.'68 x 0^'52 P.A. = 63°. 
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NGC 1333 IRAS2-A A = 2.7 mm Emission 




3''2B'°55?7 55?6 55?5 3''2B'°55f6 55f2 



Right Ascension (J2000) 

Fig. 8. — NGC 1333 IRAS2 A maps of the A = 2.7 mm contimium emission. Panel (a) is contoured 
in steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40 56.56) x the rms of panel (a) of 1.3 mJy/beam. 
Panels (b) through (d) are (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28) x a rms noise of 2.7 mJy/beam. 
(a) a = 1.3 mJy/beam; beam is 5'.'40 x 4'.'70 P.A. = 86°. (b) a = 1.2 mJy/beam; beam is 3('36 x 
3'.'16 P.A. = 45°. (c) a = 1.7 mJy/beam; beam is 1^'02 x ^'87 P.A. = 57°. (d) a = 2.7 mJy/beam; 
beam is ^.'69 x 0'.'52 P.A. = 60°. 



NGC 1333 IRAS2-B A = 2.7 mm Emission 




Fig. 9. — NGC 1333 IRAS2 B maps of the A = 2.7 mm contimium emission. Panel (a) is contoured 
in steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40 56.56) x the rms of panel (a) of 1.3 mJy/beam. 
Panels (b) through (d) are (-4 -3 -2 2 3 4 5 6 8 10 14.14) x a rms noise of 2.7 mJy/beam. (a) a = 
1.3 mJy/beam; beam is 5^'40 x 4'.'70 P.A. = 86°. (b) a = 1.2 mJy/beam; beam is 3('36 x 3^'16 P.A. 
= 45°. (c) a = 1.7 mJy/beam; beam is 1('02 x 0^'87 P.A. = 57°. (d) a = 2.7 mJy/beam; beam is 
0'.'69 X 0'.'52 P.A. = 60°. 
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SVS 13 A A = 2.7 mm Emission 



o 
o 
o 

1-3 

o 

• r— I 
-l-> 

CO 

p 




o 






O 










,100~AU 


o 

1 


ghggmgBg 


3?6 



Right Ascension (J2000) 




Fig. 10. — SVS13 A maps of the A = 2.7 mm continuum emission. All panels are contoured in 
steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40) x a rms noise of 2.2 mJy/beam. (a) o" = 1.1 
mJy/beam; beam is 5'.'40 x 4'.'64 P.A. = -70°. (b) a = 1.1 mJy/beam; beam is 3('17 x 3('05 P.A. 
= -43°. (c) a = 1.5 mJy/beam; beam is I'.'OS x I'.'OO P.A. = 57°. (d) a = 2.2 mJy/beam; beam is 
0'.'68 X 0'.'53 P.A. = 68°. The cross in panel (d) is the A = 3.6 cm peak from Rodriguez, Anglada, 
& Curiel (1997). 
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SVS 13 B A 
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Fig. 11. — SVS13 B maps of the A = 2.7 mm continuum emission. All panels are contoured in 

steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40) x a rms noise of 2.2 mJy/beam. (a) a = 1.1 
mJy/beam; beam is 5'.'40 x 4'.'64 P.A. = -70°. (b) a = 1.1 mJy/beam; beam is 3'.'17 x 3'.'05 P.A. 
= -43°. (c) a = 1.5 mJy/beam; beam is l'.'08 x I'.'OO P.A. = 57°. (d) a = 2.2 mJy/beam; beam is 
0'.'68 x 0^'53 P.A. = 68°. 
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NGC 1333 IRAS4-A A = 2.7 mm Emission 




3''39™10f6 10''4 3''S9'"10?S 10?4 



Right Ascension (J2000) 

Fig. 12. — NGC 1333 IRAS4 A maps of the A = 2.7 mm continuum emission. Panel (a) is contoured 

in steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40) x the rms of panel (a) of 3.1 mJy/beam. 
Panels (b) through (d) are (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40 56.56) x a rms noise of 2.9 
mJy/bcam. (a) o = 3.1 mJy/beam; beam is 5'.'52 x 5'.'02 P.A. = 12°. (b) a = 1.9 mJy/beam; 
beam is 3'.'02 x 2'.'81 P.A. = 1°. (c) a = 2.1 mJy/beam; beam is I'.'IS x 1'.'13 P.A. = 30°. (d) a 
= 2.9 mJy/beam; beam is 0'.'65 x 0'.'51 P.A. = 65°. The cross in panel (d) is the A = 1.3 cm peak 
from Mundy et al. (1993). 
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NGC 1333 IRAS4-B X = 2.7 mm Emission 




Right Ascension (J2000) 



Fig. 13. — NGC 1333 IRAS4 B maps of the A = 2.7 mm continuum emission. Panel (a) is contoured 

in steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40) x the rms of panel (a) of 3.1 mJy/beam. 
Panels (b) through (d) are (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40 56.56) x a rms noise of 2.9 
mJy/bcam. (a) o = 3.1 mJy/beam; beam is 5'.'52 x 5'.'02 P.A. = 12°. (b) a = 1.9 mJy/beam; 
beam is 3'.'02 x 2'.'81 P.A. = 1°. (c) a = 2.1 mJy/beam; beam is I'.'IS x 1'.'13 P.A. = 30°. (d) a 
= 2.9 mJy/beam; beam is 0'.'65 x 0'.'51 P.A. = 65°. The cross in panel (d) is the A = 1.3 cm peak 
from Mundy et al. (1993). 
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Fig. 14. — NGC 1333 IRAS4 C maps of the A = 2.7 mm continuum emission. Panel (a) is contoured 
in steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40) x the rms of panel (a) of 3.1 mJy/beam. 
Panels (b) through (d) are (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40 56.56) x a rms noise of 2.9 
mJy/beam. (a) a = 3.1 mJy/beam; beam is 5''52 x 5''02 P.A. = 12°. (b) a = 1.9 mJy/beam; 
beam is 3'.'02 x 2'.'81 P.A. = 1°. (c) a = 2.1 mJy/beam; beam is 1^'18 x 1'.'13 P.A. = 30°. (d) a = 
2.9 mJy/beam; beam is C.'eS x C.'Sl P.A. = 65°. 
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VLA 1623 A = 2.7 mm Emission 
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Fig. 15. — VLA 1623 maps of the A = 2.7 mm contimium emission. All panels are contoured in 
steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14) x a rms noise of 3.5 mJy/beam. (a) a = 3.0 mJy/beam; 
beam is 7^65 x 3'.'80 P.A. = 4°. (b) a = 2.2 mJy/beam; beam is 4'.'40 x 2'.'19 P.A. = 3°. (c) a 
= 2.0 mJy/beam; beam is l'.'44 x 0'.'74 P.A. = 10°. (d) a = 3.5 mJy/beam; beam is 0'.'95 x 0'.'39 
P.A. = 18°. The crosses indicate the A = 3.6 cm positions from Bontemps & Andre (1997). 
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Fig. 16. — IRAS 16293-2422 maps of the A = 2.7 mm continuum emission. Panel (a) is contoured 

in steps of (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40 56.56) x the rms of panel (a) of 5.8 mJy/beam. 
Panels (b) through (d) are (-4 -3 -2 2 3 4 5 6 8 10 14.14 20 28.28 40 56.56) x a rms noise of 4.8 
mJy/beam. (a) a = 5.8 mJy/beam; beam is 6'.'29 x 4'.'06 P.A. = 4°. (b) a = 4.2 mJy/bcam; beam 
is 4'.'45 x 2'.'16 P.A. = 1°. (c) a = 4.1 mJy/beam; beam is 1'.'52 x ^'76 P.A. = 7°. (d) a = 4.8 
mJy/beam; beam is 1''09 x 0^'53 P.A. = 11°. The crosses indicate the A = 2 cm positions from 
Wootten (1989). 
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100 1000 

Integrated Flux at 140 pc (mJy) 



Fig. 17. — Comparison of the ratio of the flux at 5kA and 50kA fringe spacings amplitude and the 
integrated flux of each object from Table 2, both adjusted to the distance of Taurus. The solid 
triangle symbols indicate optical/IR sources and the solid square symbols indicate the embedded 
sources. Each point is labeled with its corresponding source. 
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u,v Distance (kX) 

Fig. 18. — Model visibility curves of a 100 AU radius disk and a 5000 AU outer radius envelope 
at the distance of Taurus. These simple models assume radial power-law for both temperature and 
density (volume density, p, for the envelope model and surface density, S, for the disk model). The 
horizontal axis is antenna separation in kilo-wavelengths, and the vertical axis is arbitrary flux in 
Jy. A disk will have a 5kA/50kA ratio ~ 1, and an envelope with an outer radius of > 1000 AU 
will have a 5kA/50kA ratio > 1. 



